The interacting resonant level model in nonequilibrium: finite temperature effects 
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We study the steady-state properties as well as the relaxation dynamics of the nonequilibrium 
interacting resonant level model at finite temperatures, ft constitutes the prototype model of a corre- 
lated charge fluctuating quantum dot. The two reservoirs are held at different chemical potentials — 
the difference being the bias voltage — and different temperatures; we discuss the transport through 
as well as the occupancy of the single level dot. First, we show analytically that in the steady state 
the reservoir temperatures in competition with the other energy scales act as infrared cutoffs. This is 
rather intuitive but, depending on the parameter regime under consideration, leads to a surprisingly 
rich variety of power laws in the current as a function of the temperatures and the bias voltage with 
different interaction dependent exponents. Next we clarify how finite reservoir temperatures affect 
the dynamics. They allow to tune the interplay of the two frequencies characterizing the oscillatory 
part of the time evolution of the model at zero temperature. For the exponentially decaying part 
we disentangle the contributions of the level-lead hybridization and the temperatures to the decay 
rates. We identify a coherent-to-incoherent transition in the long time dynamics as the tempera- 
ture is raised. It occurs at an interaction dependent critical temperature. Finally, taking different 
temperatures in the reservoirs we discuss the relaxation dynamics of a temperature gradient driven 
current. 

PACS numbers: 05.10.Cc, 05.60.Gg, 72.10.Fk, 73.63.Kv 



I. INTRODUCTION 

In recent years interacting nano-structures have at- 
tracted a great deal of both experimental as well as the- 
oretical interest. They exhibit a variety of interesting 
many-body effects. Being equally fascinating and chal- 
lenging to access theoretically, these systems are subject 
of a variety of studies (for a recent review see Ref. Q]). An 
experimentally well controlled setup is given by a quan- 
tum dot with only a few electronic degrees of freedom 
contacted to a certain number of reservoirs. Due to the 
Coulomb repulsion these systems usually feature strong 
local electron correlations. A prominent consequence of 
those is the Kondo effect: if a nearly odd number of elec- 
trons reside on the quantum dot, spin fluctuations lead 
to a many-body resonance.— Besides of spin fluctuating 
quantum dot setups a typical model of interest is the 
interacting resonant level model (IRLM), which is domi- 
nated by charge fluctuations.— i 4 - The IRLM is given by a 
single level which is tunnel coupled to reservoirs (leads) 
of spinless fermions. A particle occupying the dot level 
interacts with the reservoir fermions located close to the 
dot by Coulomb repulsion. The level is coupled to an 
external gate, which allows to tune its energy. Here we 
study the IRLM in a minimal transport setup with two 
reservoirs. 

Currently the research focus shifts from the equi- 
librium to the nonequilibrium physics of quantum 
dots. One can aim at the bias voltage driven stead; 
state (which is usually assumed to exist, see Refs. 
and 0) properties,— ~— or to be even more ambitious 
at the entire relaxation dynamics from transient to 
asymptotic .— ' 10 ' 1 1 ' 1 7 ~ — In the second case it is often as- 
sumed that at some time to the quantum dot is cou- 



pled to the leads and thus relaxation from its initially 
prepared to its steady state sets in. In contrast to the 
well studied equilibrium properties, the nonequilibrium 
physics remains largely uncharted. In particular this 
holds if standard perturbative approaches fail to provide 
reliable results due to the presence of logarithmically di- 
vergent terms as it is the case in the IRLM.— For such 
problems renormalization group (RG) approache o 14 ' 21 ' 22 
might succeed to describe the system under considera- 
tion. 

Recently a functional RG (FRG) approach^ 2 - was devel- 
oped to tackle the nonequilibrium properties of interact- 
ing quantum dots 2 ^ and tested successfully in its appli- 
cation to the IRLM.— Combined with studies using 
alternative methods 2 ^— this has led to a comprehensive 
understanding of the bias voltage driven steady state as 
well as the relaxation dynamics at vanishing temperature 
of the reservoirs. By far less is known about the physics 
at elevated temperatures.— Here we aim at closing this 
gap and present a detailed study of the nonequilibrium 
IRLM with finite reservoir temperatures using FRG. We 
discuss the steady state behavior as well as the dynam- 
ics at arbitrary (possibly asymmetric) temperatures in 
the reservoirs and for a general bias voltage. Our re- 
sults are approximate, but controlled to leading order 
in the local electron-electron interaction; the RG proce- 
dure ensures that they are far superior to perturbation 
theory. Logarithms in the hopping amplitudes are re- 
summed leading to power laws with interaction depen- 
dent exponents . 3 ' 4 ' 20 ' 24 ~ 28 ' 33 ' 34 We uncover several inter- 
esting, partly surprising, temperature effects. First, we 
show analytically that in the steady state the tempera- 
tures serve as infrared cutoff scales to the renormaliza- 
tion of the tunnel couplings and as such compete with the 
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other low-energy scales [sec Eq. ([23)) ]. This is a rather 
intuitive finding, but depending on the parameter regime 
studied, the competition leads to a surprisingly rich va- 
riety of power laws in the dependence of the current 
on the temperatures and the voltage [see Eqs. ([29)) to 
(|33|) ] which can certainly not be guessed based on simple 
power-counting arguments or poor man's RG approaches. 
The interaction dependent exponents partly differ from 
the ones found at vanishing temperature. Our compre- 
hensive study of the finite temperature effects thus re- 
veals a very involved physics of the steady state current 
not reported on before. In the relaxation dynamics of the 
current and the occupancy varying the two temperatures 
allows to individually tune the amplitudes of the charac- 
teristic oscillatory terms and to vary the decay rates (see 
Fig- El) • Crucially this cannot be achieved by varying the 
two tunnel couplings to the left and right leads; we disen- 
tangle the contributions of the temperatures and tunnel 
couplings to the decay rates. As an essential intermedi- 
ate step of our analysis we present analytical expressions 
for the relaxation dynamics of the IRLM for vanishing 
interaction at arbitrary temperatures and bias voltage 
[see Eqs. (|3~i)) and ([33)) ]. To the best of our knowledge 
those were not given before. We characterize a coherent- 
to-incoherent transition as the temperature is increased. 
At this the long time dynamics changes from being expo- 
nential with overlayed oscillations to being purely expo- 
nential. This transition takes place at a critical temper- 
ature which depends in a nontrivial way on the interac- 
tion strength. Finally, we investigate the time evolution 
of a current driven by a temperature gradient instead of 
a bias voltage, which for the IRLM constitutes a so far 
unexplored setup. 

The rest of our paper is structured as follows. In Sect. 
Hll we outline the model under consideration and addi- 
tionally describe the FRG approach followed to tackle 
the problem at hand. Our result, first for the steady 
state analysis and than for the dynamics, are presented 
m Sect. Hill We conclude our paper by a short summary 
in the final Sect. [TV] 

II. MODEL AND METHOD 

A. Model 

Our model is given by the Hamiltonian 

H = H dot + Y, [H™ s + H c a oup ] . (I) 

a=L,R 

The dot part H dot consists of a linear geometry of three 
lattice sites with 

„dot TT ( n i ^ ™3\ 

H Q = en 2 -U^— + n 2 + —j 

+ (T 12 d\d 2 + T23dld 3 + ll.c.) 1 (2) 
H mt = U(riin 2 + n 2 n 3 ), (3) 



in standard second quantized notation. Here rij = d^dj is 
the occupancy operator of the spinlcss fcrmionic level j. 
The levels (sites) are connected locally through a hopping 
amplitude Ty > and a Coulomb repulsion U > 0. The 
central site can be subject to a gate voltage allowing to 
tune the onsite energy e. The second term in the single- 
particle part of the Hamiltonian is added such that e = 
corresponds to half dot filling of the central dot site 
in equilibrium. The two leads a = L,R are modeled as 
nonintcracting, 

HT = Y^A a ^- (4) 

The left (right) lead is tunnel-coupled to side 1 (3) by 
#r P = 7 Q E c L^ + H.c, (5) 

with ji = 1 and j'r = 3. For brevity we have not included 
time dependent parameters in H dot and H mt which can 
also be treated by our approach as outlined in Ref. [35| for 
reservoir temperatures T L m = 0. 

We assume that the system is prepared in a product 
density matrix state p at time to = 0, which describes a 
situation naturally arising when the coupling between the 
dot and the reservoirs vanishes for t < 0. Furthermore, 
the reservoirs are in grand canonical equilibrium with 
T a = l//3 a and chemical potentials centered around zero 
= -fJ-n = V/2 > 0, that is 

P (t=o) = p = p d ° t ® P r Q ®p%% (6) 

= e -(Hr-M/T a/Tve - ( Hr-^N a) ,T^ (7) 

where N a = Y,k a c \ °k a ■ We choose units with the Boltz- 
mann constant fee = 1, h = 1, and one electron charge 
e = I. Finally, we assume that the statistical operator 
Po ot at t = describes an initially empty quantum dot. 
This mostly studied relaxation protocol is not the one 
realized in experiments. In those the dot is initially cou- 
pled and in equilibrium with unbiased (in voltage and 
temperature) leads. At t = a bias (in voltage and/or 
temperature) is turned on.— Effects of the initial corre- 
lations present in such a setup can be studied with our 
FRG approach. This was done in Ref. [35| at vanishing 
temperature. It turned out that the initial correlations 
die out exponentially with the typical decay rate. Need- 
less to say they do not effect the steady state. We are 
confident that the finite temperature effects in the re- 
laxation dynamics discussed here equally appear in the 
presence of initial correlations. 

We aim at a limit in which the model is equivalent 
to the field theoretical IRLM .^,26-29 ,33,34 Thig ig called 

the scaling limit and can be achieved by coupling the first 
and third site much stronger to their respective reservoirs 
than to the central site. It is then that the first and 
third dot sites can effectively be incorporated into the 
lead o 2 ° l25 resulting in a single level (central site) tunnel 
coupled to and interacting with the reservoirs. We are 
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interested in the universal physics of the IRLM which is 
independent of the details of the leads band structure. 
Therefore we study structureless reservoirs (wide band 
limit) with hybridization 



T a = ivD a \j a \ 



(8) 



much larger than all other energy scales and constant 
density of states 



D a (e) = D 



-S\e\ 



(9) 



with 5 -*■ + assuring convergence of the energy integrals. 
For simplicity of depiction we will use symmetric cou- 
plings to the reservoirs of site one and three r Q = T. In 
summary the limit r^, |e|, \U\, V, \T a \ « V establishes that 
the model is equivalent to the field theoretical IRLM 
which is the prototype model of a quantum dot domi- 
nated by correlated charge fluctuations. 



B. Method 

To analyse the behavior of the IRLM at finite reservoir 
temperatures we extend the FRG methods described in 
Ref. for the nonequilibrium steady state and Ref. [25| 
for the time evolution. We here give an outline of our 
approach and discuss the essential new steps. For a more 
detailed description we refer the reader to the aforemen- 
tioned publications. 

To tackle the nonequilibrium problem at hand we em- 
ploy the Keldysh formalism.—^ In a two step proce- 
dure we first include the influence of the reservoirs on 
the noninteracting dot and secondly consider the effect 
of the two-particle interaction both in form of self-energy 
contributions. The corresponding Dyson equations are 
given as Eqs. (12), (13), (25), (29), (30), and (31) in 
Ref. HjJ The reservoir self-energy S ros is calculated 
exactly, while the contribution of the two-particle in- 
teraction is determined via the approximate FRG ap- 
proach. It was successfully applied to equilibrium 2 ? 
and nonequilibriu m 20 ! 23 ' 25 transport through correlated 
quantum dots before. In FRG one introduces a flow pa- 
rameter A in the noninteracting part of the propagation. 
We specify the cutoff procedure by coupling an auxil- 
iary structureless reservoir to each of our three dot sites 
via hybridization A. Taking the derivative of the gen- 
erating functional with respect to the A yields an exact 
infinite hierarchy of flow equations for the vertex func- 
tions which, as the only approximation of our method, 
is truncated to a given order. The cutoff-free problem is 
recovered after integrating from A = oo (where the ver- 
tices are known analytically) down to A = 0. Here we 
use the lowest truncation order; the resulting flow equa- 
tion for the interaction part of the self-energy is given in 
Eqs. (44) and (45) of Ref. I2H. The flowing self-energy ma- 
trix elements can be interpreted as flowing single-particle 
parameters. Those are the onsite energy of site one and 
three e /,A (renormalize equally in our truncation), the 



onsite energy of the central site e A as well as the hop- 
ping between sites one and two t^ 2 and between two and 
three t 23 (generically renormalize differently in nonequi- 
librium) . The problem to be solved at the end of the RG 
flow corresponds to an effective noninteracting one with 
renormailzed single-particle parameters. 

The described truncation already allows to obtain 
a comprehensive understanding of the nonequilibrium 
physics of the IRLM at zero temperature and for small 
to intermediate interactions: the logarithmically diver- 
gent terms present in lowest-order perturbation theory 
are resummed consistently leading to RG-renormalized 
hopping amplitudes featuring generic power laws with 
interaction-dependent exponents, which are correct to 
leading order in the interactioni^ - For the time evolu- 
tion FRG leads to terms exponentially decaying in time 
with interaction dependent decay rates as well as power- 
law corrections t~ K with [/-dependent exponent^ k also 
found in an alternative RG procedure . 24 ' 30 
Steady State — We employ the cutoff via auxiliary reser- 
voirs featuring infinite temperature (instead of zero tem- 
perature as in the T a = 0-study of Ref.|2Cj). We explicitly 
showed that infinite and zero temperature in the auxil- 
iary reservoirs give the same result o 25 ' 39 (the former im- 
plies technical simplifications) exemplifying the robust- 
ness of the cutoff procedure. The reservoir Keldysh self- 
energy for the present cutoff procedure takes the form (in 
lattice site space) 



//lM- 1/2 \ 

s^ cs (w) = 4;r 000 

\ f R (u) -I/2J 



(10) 



with f a being the Fermi function of reservoir a. The 
Keldysh Green function is then determined via 

G K ' A (w) = G rct ' A (w)£^ s (w)G adv ' A (w) (11) 

since the Keldysh self-energy from the two-particle inter- 
action vanishes in our truncation. We can use unchanged 
expressions for retarded and advanced Green functions as 
compared to the T = case^ the temperature does not 
enter in these quantities explicitly. The single-scale prop- 
agator S K ' A (uj) is related to the Keldysh Green function 
by Eq. (45) of Ref. [H. As the latter depends on the T a 
we have to account for finite temperature changes in the 
flow equations. Those read 



5at a = ^ f S% A (u;)du>, Tfr" = r 12 (12) 

d ^ T 23 S '23' A ( W ) da; ' T 23*°° = r 23 (13) 

^ A = -£/Ki A ^ 3 K 3 A ](^> = e (14) 

8 A ^ = JE [ S% A (lo)cLu, e' A ^° = 0. (15) 
4tt J 



Time evolution — For the time evolution we generalize 
the method outlined in Ref. [25| to finite temperatures. 
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Again we aim at the changes induced by the reservoir 
Keldysh self-energy, which enter the Keldysh Green func- 
tion G K (t,t) and the single scale propagator S K (t,t). 
For this it is necessary to evaluate integrals of the type 



/ 



dt 



smh(irT a t) 



(16) 



which can be done analytically in form of polygamma 
^(n,x) and hypergeometric functions 2-7 r i(a, b 7 c, z)4£ 
Since the formulas are rather lengthy, they are given in 
the Appendix. 



III. RESULTS 

A. Steady State 

Temperature as a cutoff scale — At this point one can pro- 
ceed with a numerical solution of the flow equations (fT2")) - 
(|T5|l . We postpone this and first report on analytical re- 
sults. To obtain those we suppress the renormalization 
of the central onsite energy e A which is 0(U 2 ) (remind 
that in our truncation we only control terms to order U) 
as well as e' A , which always appears in combination with 
the much larger scale T. This gives 



d A T2 3 



CO 

f dl[G^-\Lo)U L {u)-ll2)G^'\u)]du 

— oo 
oo 

^ / dl[G\f\^){f R ^)-ll2)Gt^)]du 

— oo 

(17) 

oo 

^ J di[Glf\co)(f^)-l/2)Gtp\u)]du 

— oo 

oo 

-EL J dl[ G f^)(f R (u;)-l/2)G^ A (cj)]du> 

— oo 

(18) 

as the remaining flow equations. We here introduced the 
star differential operator which acts only on the free 
Green function G 0,A , not on E A , in the series expansion 
qA = c ' A + G°' A J: A G - A + .. . . From these expressions we 
can extract the corresponding flow equations of the level- 
lead hybridization 9y = |t»j| /T suppressing all terms 

o(i/r 2 ) 



w 12 



UJ -e- l (A + e A 2 + e A 3 ) 

d A e A 3 =d A Q A 2 (l 3,L R). 



+ c.c. 



du> 



(19) 
(20) 



In a good approximation one can neglect the renormal- 
ization of 0i2 and 023 in the denominator — denoted as 



approximation 1 in the following — and integrate Eq. (fT9"| 
to 

dA e A 2 = -^e A 2 

xRe|^*(l,i-i^[e-/i i + i(A+e 12 + e 23 )])}, 

(21) 

with 4 r (l,x) being the trigamma function. Additionally, 
setting 1/2 + x) n 1/(2/tt 2 + x) for all Re[x] > (en- 
suring the correct value at x = as well as asymptotic 
behavior) for the right hand site of the flow equation — 
denoted as approximation 2 in the following — finally 
yields 



u 

'ttT 



9i 



2(2T L /^ + A + ei2 + e 2 3) 



' (2T l /tt + A + 0i 2 + 9 23 ) 2 + [(e - Ml)/2] 2 ' 

(22) 

After reintroducing the ultraviolet cutoff 4 ^ T one can in- 
tegrate this equation analytically and obtains 



B A=0 



-2UI(ttT)+0{U 2 ) 



\e-fi L \,T L «T K «T 



V -2UI(*T) + 0{U>) \ e \ : T K ,T L «V«T 
| e |-2t//(-r)+0(i/ 2 ) V t k t, « | e | « r 



-2UI(kT)+0(U 2 ) 



K,J-L 

e-[i L \,T K «T L «F 



(23) 



with r = T12 + T23 and Tk being the universal equilibrium 
low energy scale of the modet^ ^ 29 ' 33 ' 34 defined via the 
charge susceptibility- 2 ^ 



dn 
X= lTe 



c=0 



ttT k 



(24) 



Equations analogous to ([2"TT) - ([2"3"|) arc obtained for the 
renormalization of 023° by replacing (1 -> 3, L -*■ R). It 
is apparent [see Eq. (|23|) ] that the reservoir temperatures 
provide infrared cutoffs. Equation (|22j) indicates the in- 
terplay of the different cutoff scales e - /j, a , ©12 + @23,T a 
entering in the denominator explicitly. We find that the 
left (right) hybridization is cut off by the temperature 
of the left (right) reservoir only. This invokes important 
consequences for the behavior of the current outlined in 
more detail in the next section. Different cutoff scales for 
the left and right hybridization have been observed be- 
fore (at T a = 0); it is |e-/ij,| (|e-/ifl|) and not V = (J>r 
which enters as a cutoff for 812 (©23)^ 

We next elaborate on the quality of the approximations 
1 and 2. To this end we take the full numerical solution 
of the flow equations (the only approximation being the 
truncation) and step by step apply approximations 1 and 
2. A comparison for the renormalized line width A=O as 
a function of temperature is shown in Fig. [TJ For the 
entire temperature range (Tj, = Tr = T and T12 = T2 3 = r 
for simplicity) the results of approximations 1 and 2 agree 
well with the numerical solution. In particular for T 
Tk, the power law Eq. ([23| with the correct exponent is 
reproduced on each level of approximation. 
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FIG. 1: (Color online) Renormalized line width A as a 
function of temperature T. Parameters are chosen as e/Tfc = 
V/T K = 0, T12/T = r 23 /r = 0.0025, T = T L =T R and U/T = 0.1. 
Since left and right lead are chosen symmetrically Q12 = O23 = 
<d holds. Notice the double logarithmic scale, which indicates 
the power law behavior for large T. The inset shows the 
replacement of the trigamma function used in approximation 
2. 



Observables — Within our truncation the effect of the in- 
teraction can be cast into effective single-particle param- 
eters. To gain analytical insights it is therefore sufficient 
to derive expressions for the occupancy and the current 
at vanishing interaction. One can simply supplement the 
noninteracting formulas with the effective parameters ob- 
tained at the end of the flow to incorporate the interac- 
tion. We focus on the current which is the most interest- 
ing observable in our transport setup. We start from the 
Meir-Wingreen form 4 ^ 



rct,A=0|2 
13 



J L =4r 2 J du[f L {u)-f R {w)]\G 



(25) 



where terms 0(U 2 ) (which we do not controll in any 
case) as well as 0(1/T 2 ) (scaling limit) were neglected in 
the second step. The integral can be solved and yields 



Jl=4 



^12 ^23 



9f 



A=0 + A=O 



23 a=L,R 



X Im 



*(i-^W^- e )-(e^° + e^ ))) 



(26) 



with ty(x) being the digamma function, s(L) = -1 and 
s(R) = 1. Combining this expression with the flow Eq. 
(|2"l"|) and its analogue for ©2 3 allows us to obtain a com- 
prehensive analytical understanding of the steady-state 
current. 



In the ohmic regime (in which the current depends 
linearly on the voltage), for Tl = Tr = T, t 12 = r 2 3 = r, 
and at low T where Tjc » T, V,\e\ the T-dependence of 
the current was earlier studied in Ref. using a field 
theoretical approach. In this limit we find 



e 



A=0 



e 



A=0 



cxp 



2U 


if 2 


7T 4 y 


_~rv 


{ 6 


60 / 



(27) 

where we neglected terms 0(T e ) with T = nT/T K . Plug- 
ging the renormalized hybridizations into the current for- 
mula (|2"6"|) yields 



■lL = l-lf 2 
V 3 



(7_ _ 2U \ 
1 15 9T7T/ 



T 



(28) 



up to 0(T 6 ) and 0{U 2 )._ The ratio R = g 4 /g% of the 
prefactors of T 2 (<? 2 ) and T 4 (54) is given by R= 21/5- 
2U/Tir, which is the same result as found in Ref. [2fj| (see 
also Ref. S3). 

At T L j R = 0, for left-right symmetric tunnel couplings, 



e = and for sufficiently large bias voltages Eq. (|23j) 
leads to a power-law suppression Jl ~ V~ v of the cur- 
rent with v = -^^2S It was later shown^ 4 ^ that re- 
laxing the second and/or third requirement leads to a 
more complicated form of the current-voltage character- 
istics. Allowing for a left-right asymmetry in the reser- 
voir's temperatures introduces another level of complex- 
ity. We break it up by considering the symmetric case 
T a = T » \e - fj, a \,T}(, first. In this regime the flow of 
both level-lead couplings is cut off by the temperature 
and the current reads 



(29) 



as long as V + or e + 0. We emphasize that the expo- 
nent is different from the above given v appearing in the 
the voltage dependence. To the best of our knowlegde it 
has not been found before. For V = and e = 0, Jl = 
independent of the temperature difference in the reser- 
voirs. This follows right away from half filling. In the 
regime T L » Tr » \e - pL a \, T K , 



JL~(T L T a )-*-*- C2TL + TR 



(30) 



where ci stems from the asymmetry in the renormal- 
ization of the left and right hoppings (the exact form 
of c\ is irrelevant for the following discussion) and c 2 = 
-(e-V/2)/(e + V/2). Off resonance (e * V/ 2) the current 
is given by 



Jl - (T l TrY 



t: 



L + ^ T R T 

but on resonanc o 24l3 ° (e = V/2) it changes to 

Tr 



Jl - (TlTrY 



T, 



+ CiT ; 



(31) 



(32) 



R 
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Similarly one finds for T L » V » T R , \e\,T K , 

j l Ku L — • ( 33 ) 

This exemplifies how the temperature significantly 
changes the qualitative behavior of the current, although 
it enters the renormalization of the hybridizations in an 
("simple") intuitive way (as an infrared cutoff). We em- 
phasize that the subtle interplay of the lead tempera- 
tures and the voltage revealed by the above equations 
can only be uncovered using a method (such as ours) 
which allows for an unbiased RG-likc treatment of mul- 
tiple energy scales. 

B. Time evolution 

The noninteracting case — To gain a comprehensive un- 
derstanding of the effect of finite reservoir temperatures 
on the relaxation dynamics of the IRLM it is advanta- 
geous to first study the noninteracting case. We derive 
closed analytical expressions for the time dependence of 
the current and occupancy, which to the best of our 



knowledge were not presented before. Based on those 
the effect of T a > on the relaxation rates and charac- 
teristic oscillations can be worked out in detail. A com- 
parison with the U > results obtained numerically by 
solving the FRG flow equations and computing n(t) and 
Ji,(i) then allows to assess the correlation effects. The 
correlation physics is particular prominent and transpar- 
ent for |e ± V/2\ » Tk- It was earlier shown that for 
T a = in this regime the relaxation dynamics for U > 
and sufficiently large times can be obtained by replacing 
the time-independent bare single-particle parameters in 
analytical U = expressions by the time-averaged renor- 
malized ones.— We show that the same holds for finite 
reservoir temperatures. This constitutes a second reason 
for the U = expressions of n(t) and Jl(^) providing the 
basis of our understanding also of the interacting model. 

We introduce the hybridization of the single-level 
model Tij. To connect to the three site model in the 
scaling limit one has to set = |r i:( | 2 /r. For simplicity 
we focus on the case of symmetric couplings to the reser- 
voirs = r. The occupancy n{t) with the initial value 
no reads (for expressions of the occupancy in terms of 
Green functions see Ref. [25T ) 



i(f) 4"^ 4f ' (i " 2fi ° )+ ^ r ° im i 



(i(e-/j Q ) -2f -7rT< 



2ttT„ 



/ nT a \ 2nT a J 



2 e (i(e- i "a)-2f-7rT a )t 

(i(e-/i a ) - 2f--nT a ) 

2 e We-Ma)-2f-7rX a )t 

(i(e - fi a ) + 2T -TiT a ) 



2-^1 



-(i(e-/x a )-2f-7rr a ) -(i(e - fj, a ) - 2T - nT a ) _ 2vTat 

J-i i +1,6 



27rT Q 



2ttT q 



-(i(e-/j, a ) + 2T -7rT a ) -(i(e-fj, a ) + 2T -7rT a ) - 2 7rT„* 
1, , + l,e 



2ttT„ 



2ttT„ 



(34) 



r 



Interestingly, the rate 4r appearing in the first exponcn- tions expressing the current in terms of Green functions 
tial term remains independent of the temperatures. For we again refer to Ref. |25| ) 
the time dependence of the current we find (for equa- 



m -1 - 2n(t) + 4T^ ^ -1 

r 2ttT l 

AT L Re\i = 2^1 

i(HL -e)+2T + ttT l 



i(fi L - e) + 2T + nT L -(i(e - fj, a ) + 2T - nT a ) - 27r T„t 
1, — = , ^rz + l,e 



2ttT q 



2irT a 



(35) 



Up to the above mentioned first time dependent term 
in n(t), which remains unchanged if T a is increased, 
the dynamics is affected by the reservoir temperatures 



[e.g. compare to Eqs. (109)-(111) of Ref. HJ. The 
long-time behavior of the remaining terms in n(t) and 
J hit) changes from being governed by an exponential re- 
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laxation in combination with a power-law correction for 
T = 24 ] 25 i 30 to an infinite series of exponential terms 
with temperature dependent rates 2T + 2nirT a and n e N 
at T > 0. The frequencies characterising the oscillatory 
part of the behavior of fi(t) and Jh{t) are the same as for 
T = O i 25 ' 30 in general one finds the frequencies |e ± V/2|, 
where the amplitude of |e + V /2\ (|e - V/2|) is suppressed 
in the regime |e ± V /2| » Tk in the current Jl (Jr)- 

Obviously the reservoir temperatures do not play the 
same role as the hybridisation (remind that Ti2 = 1^23 = 
T). In contrast to the latter, which enters in every relax- 
ation rate of Eqs. and the temperature allows 
to tune the influence of the different terms with respect 
to each other; the rates dependent solely on either Tl or 
on Tr. Furthermore, even with asymmetric couplings to 
the left and right reservoir, that is generalizing Eqs. f34|) 
and ([55)1 . only the sum of the hybridisations enters the 
decay rates. The temperatures instead allow to tune the 
relaxation rates imprinted by the left and right reservoir 
independently. 

The interacting case — In the developed FRG approach 
we determine the Keldysh and retarded Green functions 
of the interacting system. Given those we can numer- 
ically compute the observables as outlined in Ref. I25I . 
We postpone this and first consider the time dependent 
rcnormalization of the single-particle parameters. When 
combined with Eqs. and the latter give valuable 
insights in the effect of the two-particle interaction. We 
carefully checked that the numerical parameters underly- 
ing the (numerical) solution of the flow equations (for de- 
tails see Ref. HH) are chosen such that the results are com- 
pletely converged on the scales of all figures shown. In 
fact, it is one of the advantages of our approach (partic- 
ular, compared to purely numerical ones) that for small 
to intermediate U higly accurate results can be obtained 
for arbitrarily large times. 

For arbitrary T a and the relevant times T^t g; t/T 
the renormalization of the onsite energy e A=0 (i) - e = 

^22 t A ~°(^) i s °f or der U 2 and can safely be neglected. 
The time dependence of the renormalized hopping am- 
plitudes is shown in Fig. [5] (for T12 = T23 = r). It depends 
only weakly on the temperature for the considered case 
of large |e± V /2\. As for T = 0^ the effective hopping am- 
plitudes quickly (on a scale T^-i ~ t/T) start to oscillate 
around their steady state values with frequencies |e— V/2| 
and |e+V/2| for r-^f and t^ respectively. Note that due 
to the left-right asymmetry induced by the voltage and 
different reservoir temperatures the renormalized level- 
lead couplings gcnerically differ even for the bare ones 
being equal. Increasing the temperature of the left (right) 
reservoir suppresses oscillations in the renormalized r-^f 
( T 23 = °)i this will manifest as an interesting tunability in 
observables (see Fig. [5]). For the sake of simplicity we 
only show the real part at this point, but the same holds 
for the much smaller imaginary one as well. 

The analysis of the time dependence of the renormal- 
ized single-particle parameters in the limit |e± V/2\ » Tk 
shows that to interpret the results obtained by the nu- 



Re[\~ -T]/r K : T L /T K =0.5; T R /T K =0.5 • Re[T, ° -l]/T K : TJT^OS; T R /T K — 0.5 

RelT^-T]/^: T L /T K =0.7; T R /T K =0.3 ReK^-Xl/T^ T L /T K =0.7; T R /T K =0.3 

Refu, 2"°-lFr K : T L /T K =0.9; T R /T K =0. 1 Re[x*"°-i:]/T 1( : T L /T R =0.9; T R /T K =0. 1 

~r\^_~ ' T" 
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■ . >.;•■:■■.'■:' •:■ '■>:■ ■:.<-■:.:".>:. ■ »i ■•■ ■•■ ■:•> ■•> ■•■ >; ->s .m/.Mj'S 




FIG. 2: (Color online) Time dependence of the renormalized 
hoppings 7"i2 =0 - t = E^2 t,A ~° £ C between sites 1 and 2, and 
T^f - r = Ejg*' ~° e C between sites 2 and 3. The parameters 
are: t 12 = r 23 = r with t/T = 0.0025, e/T K = V/T K = 10, and 
U /T = 0.2. The insets show a zoom into the regime where one 
can most clearly distinguish the temperature's effect on the 
oscillations' amplitude. 



merical solution of the FRG Sow equations and numer- 
ical computation of the observables for Tjct ^ t/T one 
can refer to Eqs. (|34j) and (|35|) with the bare parameters 
replaced by the time-averaged renormalized ones. More 
precisely we have to extend these expressions to the case 
of left-right asymmetric hybridzations. For briefness we 
do not give those here but note that the ratio of the hy- 
bridizations enters as prefactors in the different terms. 

We now discuss our numerical FRG results for n(t) 
and Jh{t) at U > 0. To understand those we use the just 
established relation to the noninteracting dynamics. Fig- 
ures [3] and 2] show the time evolution of the occupancy 
and the current in the limit |e ± V/2\ » Tk for different 
U . Increasing the temperature drastically suppresses the 
amplitude of the oscillatory terms. The quality factor 
is decreased by the increased decay rate. One can tune 
the quality factor of the two frequencies |e ± V/2\ inde- 
pendently of each other, which is not possible via the 
hybridization. Figure [5] illustrates this point for the time 
dependence of the occupancy. For vanishing tempera- 
ture gradient AT = it shows a superposition of the two 
frequencies |e ± V/2\ (bottom curve in Fig. [5]). Raising 
the temperature gradient pronounces the contribution of 
the frequency |e + V/2\ belonging to the colder reservoir 
until the signal appears almost sinusoidal (ovcrlayed by 
an exponential decay) at maximum temperature gradi- 
ent (top curve in Fig. [5]). The inset shows the imaginary 
part of the numerical Laplace transform (normalized to 
its largest value) of the occupany for the smallest and 
largest temperature gradient. The positions of the fre- 
quencies |e ± V /2\ are indicated by the vertical arrows. 
With increasing AT the feature at \e-V/2\ is suppressed, 
while the one at |e + V/2\ is enhanced. 

Next we study the far-from-equilibrium case V » 
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FIG. 3: (Color online) FRG data for the time evolution of 
the central-site occupancy. The parameters are the same as 
in Fig. [2] and symmetric temperatures T a = T. The arrows 
right to the graph indicate the steady-state values obtained 
by the nonequilibrium steady-state functional RG discussed 
earlier. 
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AT7T K =0.2 

AT7T K =0.4 

AT7T K =0.6 

AT7T K =0.8 

AT7T =1 



increasing AT 
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FIG. 5: (Color online) The same as Fig. El with U/T = 0.2, but 
featuring a temperature gradient Tl = T + AT/2 and Tr = T — 
AT/2 for T = 0.5Tk- We find how the temperature gradient 
pronounces the frequency |e + V/2\ belonging to the colder 
reservoir. This can also be found in the inset, which shows 
the numerical Laplace transform with the positions of the 
frequencies \e±V/2\ indicated by arrows. 



0.02 



U=0; T7T =0 




FIG. 4: (Color online) The same as in Fig. El but for the 
current leaving the left reservoir. 



Tx,e = 0, which on general grounds is the most intrigu- 
ing one. More precisely we consider the limit V » T^,T 
with Tl = Tr = T and e = 0. At e = the expressions for 
the time evolution are particularly simple and for suffi- 
ciently large times 7rTt » 1 and small interactions we 
obtain 



Jh(t) = Jstat + e 



-46" 



2TIm 



4e (-iV/2-2e A - -7rT)i 

iV + 49 A=0 + 2ttT 



with the renormalized steady-state values 



A=0 



(36) 
A=O - 

O^f u and the stationary current J s tat- We note that for 
initially equal 0y and e = also the renormalized hy- 



bridizations remain equal even for V + 0. From Eq. (|36[) 
we can read off a coherent-to-incoherent transition in the 
long-time dynamics at temperature 

29 A=0 . 



irT c 



(37) 



The long time dynamics switches from being exponential 
with an overlayed oscillation (second term of Eq. ([3"6"]l 



being the dominate one at large times) for T < T c to 
purely (monotonic) exponential relaxation (first term of 
Eq. ([3^1) being the dominate one at large times) for T > 
T c . To determine T c we can take the expression for A=O 
corresponding to the considered parameter regime V » 



T K ,T of Eq. 
obtain 



T 



A=0 



f(v) 



2U/7r/r 



2r z / T \ 



2U/n/r 



With this we 



(38) 



for the critical temperature at which the crossover from 
coherent to incoherent behavior occurs. Relaxing the 
condition of equal temperatures in the left and right 
reservoirs but keeping V » Tk,Tl,Th the transition is 
found at min(Tj / ,TR) = T c with the same T c as defined 
above. 

We close the discussion by considering the time de- 
pendence of a current which is not induced by a finite 
bias voltage but by a finite temperature gradient accross 
the quantum dot. To obtain a nonvanishing steady-state 
value of the current one needs to choose e + 0. Figure [6] 
shows the current for this parameter regime. Note that 
the only frequency in the time evolution is given by e 
as V = 0. The interaction strength enhances the ampli- 
tude of the oscillations, but decreases the steady-state 
current (if measured with respect to Tk)- This inves- 
tigation opens the road to thermal and thermoelectric 
transport studies which will be presented elsewhere. 



IV. SUMMARY 

We extended the functional RG analysis of the 
nonequilibrium interacting resonant level model in the 
scaling limit, which is the prototype of a simple charge 
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FIG. 6: (Color online) Time evolution for the temperature 
gradient induced current leaving the left reservoir J Lit). The 
parameters read t/T = 0.0025, e/T K = 10 and V/T K = 0. We 
choose Tl = Tk and Tr = (choosing < Tr « Tl does not 
alter the qualitative behavior). 

fluctuating quantum dot, to the case of arbitrary reser- 
voir temperatures. Our study includes all time regimes 
from transient to asymptotic with the ultimate limit of 
the steady state. All results developed are controlled 
to leading order in the interaction, but known to be far 
superior to simple perturbation theory as our RG pro- 
cedure implies the resummation of logarithmically diver- 
gent terms to power laws. 

First we discussed the role of temperature as an RG 
cutoff competing with the other energy scales of the 
model in the steady state. Although the temperature 
enters the renormalization of the hybridization in an intu- 
itive way, the current was shown to exhibit a vast variety 
of different power laws (in temperature and bias voltage) 
with interaction dependent exponents for different pa- 
rameter regimes. We than clarified the role of finite tem- 
peratures for the relaxation dynamics. We focused on the 
most transparent case |e± V/2| » Tjc,Tl, Tr in which the 
physics of the IRLM can be understood by interpreting 



the interacting system as an effective free one featuring 
renormalized parameters. The finite reservoir tempera- 
tures T a (only) partly act similar to the hybridizations; 
increasing T a increases the relaxation rate. Our detailed 
analysis shows that this obvious effect does not exhaust 
the role of temperature in the relaxation dynamics. We 
have shown explicitly that one can use the T a to suppress 
the oscillatory contributions to the occupancy with the 
frequencies |e ± V/2| individually which is impossible by 
varying the hybridizations. We charcterized a coherent- 
to-incoherent transition in the long-time relaxation dy- 
namics as the temperature is increased. The critical 
temperature at which the transition occurs depends in 
a nontrivial way on the two-particle interaction. Finally, 
we discussed a current, which is not driven by a bias volt- 
age, but by a temperature gradient. We emphasize that, 
similar to the T = case, our approach can also be used 
to tackle explicitly time-dependent Hamiltonians, which 
constitutes a field of current interest. A direction of fu- 
ture research which opens up directly from the present 
work is the field of thermal and thermoelectric transport. 
In this context quantum dots are promising candidates 
for highly efficient energy conversion devices and thus 
have recently gained sizeable interest.— ~— 
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Appendix 

In the main part of this work it was discussed, that one 
needs to calculate the Keldysh Green function G K (t,t). 
It is determined by the integral 



J d Sl J d S2 G rct (t„ +1 , Sl )I]K s ( Sl ,,s 2 )G adv (.s 2 ,t m+1 ) (A.l) 



with S res (si,s 2 ) given by 



££(t',t)=r 



/-Tie - *"* (t ' _t) 



sinh[TTT L (t'-t±iS)] 1 



-T R e 



-i[i R (t'-t) ■ 






-. 1 

•>± sinh[7rT R (i'-t±i5)] / 



(A.2) 



The retarded Green function can be evaluated as in the with Res^,; and u>i being the residues and the poles of 
T = case via 

Gf(t,t') = -z^Res^e-^"'), (A.3) u - (h** + E «t _ iA + Ep t) " 

i=i 



10 



11 



i.i 

22 12 



13 



1 1 + 



n 2 



|.A ,2 







(ui 1 -ui 2 )(ui 1 -ui 3 ) 



^2-^1 T 12 



r 12~23 



(w2-wi)(iJ2-aJ3) w 2 -w,s W2^3 (^2 -"1 ) (w 2 -UJ 3 ) 



"3~"1 T 12 



33 



Resii, n (r^ -» t 2 a ) 



21 ReBia,nWa - WU*) 

31 Resi3,n(ri A 2 ,r^ -> (t^)*, (t^)*) 



23 



Resi2,n(r 12 -» r 23 ) 



TABLE I: Residues of Eq. (fA~4)) 



32 Res 23 ,„(r 2 A 3 ■+ (r 2 A 3 )*) 



Introducing the time dependent effective parameters 



T 23 = T + ^^23 a U° ws to express the poles as 



^ = e' A - i(T + A) , lu 2/3 = i (e A + e' A - ir - 2iA T \/-(r - i£ A + ie'> A ) 2 + 4|t£| 2 + 4|t 2 a 3 | 2 j . (A.5) 

The corresponding residues are given explicitly in Tab. [TJ 

For the integral Eq. (|A.1[) one substitutes the 'center of time' T = t\ + t 2 and the 'relative time' At = t%-ti: 



efts 



im+1-ti.+i 2i„ +1 +Ai t m -t n 2t m+1 -At t m+ i-t„ 2t m+ i-At 

J dAt J dT+ J dAt J dT+ J dAt J dT 



2t m -At t m+1 -t„ 



2t m -At 



2t„+At 



(A.6) 



for t n+ i - t n > t m+ i - t m . The opposite case follows analogously. It proves advantageous to separate the problem into 
n = m and n + m. We start with the case n = m. With the above substitution one can write the i,j matrix element 
of Eq. (|A.1[) in terms of digamma ^/(x) and hypcrgometric functions 2 .Fi(a, b, c, z)£L 



•i + l Era+l 



J dsi J ds 2 G 10t (t n+1 ,s 1 )Y,f ca (s 1 ,s 2 )G adv (s 2 ,t n+1 ) 

Lt„ t„ 

-S 2t„ + i+At t„+i-t„ 2i„ + i-Ai 

J dAt J dT+ J dAt J dT 

t„-t„ + i 2t n -At 5 2t n +At 



lim V TaReSjQ fc Res* a ; e 
<s-»o 



-iAtii fci t„ + i 



sinh(7rT Q Ai) 



-L,R 
k, 1=1,2, 3 



gijLtaAigiiTAufeig-iiAtCwfe+wf) 



fc ,1=1,2,3 



2ttT q I V 



-ifi a + iujk + ttTq \ / - + irT a 



2irT n 



) \ 2nT a ) 



D H Ma +fa;»-MrT a )(t»-Wi) / -i^+icjfc + ^Ta -l^ a + lCJ k + TlT a 2*T a (t n -t n + 1 )\ 
2-^ 1 I 1) ; r 1,6 v y I 

+ iwfe + 7rT Q \ 2?rT a 2nT a } 

i(i a - iujf + nT a i(i a - ilui + nT a i 1 ^T a (t n -t n+1 ) 



-ifj, 



^(•ijLia— iw ! *+7rT CI )(t n -i„ + i) 



i/io, - + 7tTo, 



JAui k it n 



1 



2^l(l,- 

+ iu; ; * + 7rT Q \ / - iw fc + -nT a 



2nT a 



2nT a 



') 



27rT Q l V 2ixT a 

-2 



/ I 2nT a ) 



-ilia + icjf + ttT c 



H ' 27rT a ' 2irT a ' 6 J 



(i ( u a ,-iw l! +7rT a )(t n -t„ + i) 
+ - = %?1 



7T f-[ ~ iuJ k + »Ma ~ V£k + T^rq | - | fi 27rT„(t„-t„ + i) ' 



27rT Q 



2ttT q 



)])• 



(A.7) 
(A.8) 



where one has exploited that for b > 



lim 



7-2Tl 1, 

a + 6 V 



26 ' 26 



c + 6 



26 26 
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and defined 



Au)M = UJk - U)f 



In the indices of the residues in Eq. (|A.8|) one has to replace a = L by 1 and a = R by 3. For the case m * 
analogously finds 

p e i/i„At e iiTAw fc , 



(A.10) 
n one 



y dsi y ds 2 G ret (t„ +1 , sl )sf es ( S1 , S2 )G adv ( S2 ,t m+1 ) 



1 

ij k, 1-1,2, 3 



y ^ y 



tm+l-Wl 2t„ + 1 +At t m -t„ 2t TO+ i-At tm+l-tn 2f m+ i-At 

y y dT 

t m -*n+l 2t m -At tm+l-Wl 

4T r 

X T a Res ia , k Re S * a 



y dAt y 

2t m ~-At t m -t 



2t„+At 



sinh(7rT Q Ai) 



|_ e -jiAt(w fc +aj,*) 



M-1,2,3 



p(VQ-'i"r+ 7r ' r o)(*m + l-tr 1 + l) / 7 ',, _ ,V)» + —T' _ J/J* + 7rT 

_ iAu> k it n+1 V Til 1 a "° * ' 

ljJL a - IUJ? + TTl a \ 



2nT a 



■ iwl + irT a e 27rT a (i m+1 -t„ +1 )\ 

2nT a ' j 

in a ~ iu* + TrT a \ ' 27rT Q 27rT Q / 

l >e 2»T-(*»-Wi)J 



(i^ a --»ujfc+'7rT a )(t m -t n+ i) 



(l,^ 



2ttT q 
- iw fe + 7tT q 



2wT a 



2TrT n 



l^a - ZWfc + 7rJ Q 

A iA, M ^ 1 e (j "°^ wr ° XW " t ° +l) ~ /, - ^fc + ^ */i a - ico k + nT a 

+ e + 2-0 1 1, j + 1, 

7rT Q \ 2ixT a 2irT a 



i[i a - iuj k 



Ji^ a -iui k ^T a )(t m -t n ) 
+ iAu> kl t m c 



iAukltm+i £ 



e 27rT Q (t m+1 -t„ +1 )\ 

+^T„)(t m -t„) / i/z Q - jgjfc + 7rr a ifj, a -ioj k + irT a | ^ e 27rT a (t m -t„)\ 
ia;^ + 7rT a ~ \ ' 2irT a 2TiT a j 

^T a )(t m+1 -t„) ^ ^ i^ a -iuj k + 1 TT a in a - ioj k +TtT a | 1 2 „ Ta (tm+1 -t n ) \ 

2nT a + ,B j 



_ / ifx a - iuj k + nT a i^ a -iuj k +TrT t 
; — — 2-ri I 1, 

lfl a - lUJk + 7T 1 a \ 



2nT a 



- 1 - a v ** *' en ^ " Q: 

i/i a - iuif + irT a " \ ' 2i^T a 2i:T a J 

+ 7rT Q \ ' 2iiT a 2iiT a j 



ifj, a - iu 



2nT a 



Similarly one can derive S K (t,t), which additionally involve s the trigamma function ^ ( 1 , : 
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